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We show that randomly choosing the matrices in a completely positive map from the unitary 
group gives a quantum expander. We consider Hermitian and non-Hermitian cases, and we provide 
asymptotically tight bounds in the Hermitian case on the typical value of the second largest eigen- 
value. The key idea is the use of Schwinger-Dyson equations from lattice gauge theory to efficiently 
compute averages over the unitary group. 
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Recently, twopapersP, |3| introduced the idea of expander maps: quantum analogues of expander graphs. An 
expander graph[j| may be defined in several ways. One is the property of having a large number of vertices, a 
■ small coordination number for each vertex, and also having a gap in the spectrum of the diffusion equation on the 
[ graph, so that a particle classically diffusing on an expander graph rapidly loses memory of where it started. In the 
quantum case, we replace the random process of diffusion by a completely positive, trace preserving map £{M). We 
define a quantum expander to be such a map from the space of N-hy-N matrices M to the same space with the 
following properties. First, N is large, in analogy to the large number of vertices. Second, the map has eigenvalues 
Ai, A2, A3, A^v^ , with Ai = 1 and |Aa| < 1 — S for all a > 1 so that the eigenvalue spectrum has a gap. Finally, the 



map can be written as 
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Ph; £{M) ^Y.A\s)MA{s) (1) 

c ; 

^ , for some relatively small value of D, with = 1 so that the map is trace preserving, and with 

X^fLi A^{s)A{s) = 1, so that the eigenvector corresponding to eigenvalue Ai is {l/VN)t. Here, 1 is the iV-by-iV unit 
I— ' matrix. This requirement of small D is in analogy to the low coordination number. 

These maps were applied in [l| to construct many-body states in one dimension with the property of having a short 
correlation length (this corresponds to the gap S in the spectrum of eigenvalues of £), small Hilbert space dimension on 
' each site (this corresponds to the small D) , and yet large entanglement entropy (this corresponds to the large entropy 
I of the eigenvector of £ with unit eigenvalue) . Since expander graphs have a large number of applications in problems 
' dealing with classical statistics, such as in error-correcting codesQ, derandomization, and the PCP theoremfG], to 
name a few, it seems worth further exploring the quantum case. 
\^ . A number of possible forms of an expander map are possible: in [2] an expander was defined as having 
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for some unitary matrices U{s), and in fact this is the form of A(s) considered in this paper. However, the more 
general definition with arbitrary A{s) constrained by ^fii A^{s)A{s) = l,Er=i A{s)A^s) = 1 seems also useful; in 
fact, although we do not consider it in this paper, it may be useful to weaken this constraint further, and explore the 
properties of completely positive, trace preserving maps, with no other constraint on the A, requiring only that the 
entropy of the density matrix which is the eigenvector with unit eigenvalue is large Q. 

Our goal is to try to find families of maps with arbitrarily large N, such that the gap 5 is bounded below by some 
A'^-independent constant and such that D does not grow too rapidly with N. The first paper[lj| provided an explicit 
construction of such a family of maps with D of order log(A'^) and provided numerical evidence for an alternate 
construction with D independent of A^. The second paperQ gave yet a different construction with D of order log(A^) 
but also provided a construction that had D independent of A^ and succeeded in proving an A^-independent lower 
bound on the gap S in this case. 

Experience with expander graphs suggests that, while finding deterministic constructions of them is difficult 0] , with 
high probability a random graph of fixed coordination number greater than 2 is an expander^]. Thus, the natural 
question is to investigate whether Eqs. (jll2p will give an expander map if the matrices U are chosen randomly from the 
unitary group U {N) using the Haar measure. We consider two cases. In the first case, the map £ is non-Hermitian and 
the D matrices are chosen independently at random. In the second case the matrices U (s) are chosen independently 
at random for s — 1...D/2 and we pick U{s + D/2) = U{s)^. In this case, D is even, and the map £ is Hermitian 
and has real eigenvalues. In this paper we begin in generality with the non-Hermitian case, but then restrict to the 
Hermitian case for simplicity of notation. 
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In the Hermitian case, we consider D > 4, while in the non-Hermitian case we consider D > 2, as otherwise we 
would clearly not have an expander. Let A2 be the eigenvalue with the second largest absolute value of all eigenvalues 
other than Ai. Let 



A» ^ (3) 

The main result of this paper is that, in the Hermitian case, for any e > the probability that [A2I is within e of 
Xh approaches unity as ^ cxd. Interestingly, this is the same as the recently proven tight bound [9] in the classical 
case, but the proof in the quantum case is much simpler. 

The proof is based on a version of the trace method. We begin by introducing the trace method and describing 
its application to the Hermitian and non-Hermitian cases. We then give lower bounds on IA2I based on the return 
probability of a random walk on a Cayley tree and discuss some numerical results. We next introduce a set of 
Schwinger-Dyson equations, analogous to those used in lattice gauge theory This is the key step which enables us 
to take averages over the unitary group efficiently. We will use these equations to develop a convergent perturbation 
theory in 1/iV for various traces of unitary matrices, and bound the correction terms in this perturbation theory. We 
start with a loose bound, giving a loose bound on IA2I, and then tighten to get the tight bound above. Finally, in an 
appendix we discuss a related problem of "quantum edge expanders" , which gives an analogue in the quantum case 
of the combinatorial definition of a classical expander graph. 

The space of N-hy-N complex matrices M has a natural inner product: iM,N) = tr(M'fiV). With respect to this 
inner product, an orthonormal basis of matrices consists of the matrices M{i,j), defined to have a 1 in the i-th row 
and j-th column, and zeroes everywhere else. Given this inner product, we can consider the space of N-hy-N matrices 
as an iV^-dimensional vector space, with £ acting as a linear operator on this space. Then, in the Hermitian case, it 
is possible to find a linear operator V, which is unitary with respect to this inner product, such that £ = V^AV, with 
A being a diagonal matrix with entries Xa- Note that here £,V, and A are all N'^-hy-N'^ dimensional matrices. In 
the non-Hermitian case, we can write £ = V^TV, with T an upper triangular matrix whose diagonal entries are the 
eigenvalues Xa- Thus, 

^(£™(M(^,J)),£™(M(^,J))) ^ J] (T"(A^(^, j), T™(M(z, j))) (4) 

2m 
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a=l 



where £"^{M) denotes acting with the map £ successively m times on M, and similarly for T"''{M). In the case where 
£ is Hermitian, Eq. ^ is an equality. 

To simplify notation, we now restrict to the Hermitian case. In this case, Eq. ^ can be replaced by 

^(m(», j),5™(A/(*, j))) = i^-i" ^ 1 + i^^r, (5) 

i,j a=l 

where we pick m to be an even integer. Then, 

a=l i,j 

{j^Y H- E ^<U(Srn+ D/2)...U{S2 + D/2)U{si + D/2))iY{U{si)U{s2)...U{s,n))]. 



For notational convenience, we identify Sj 4- D with s.^ throughout this paper, so that Si is a periodic variable with 
period D. 

Averaging U{1), ...,U{D) over the unitary group we find that i;^^^^. (M(i, "(M(i, j)))] = i^E^^i l^al™], 
where E[...\ denotes the given average. Averaging Eq. ^ we find 

^^^(l^Til E- E E,{s,,...,s,^)^E[Y\Xari (7) 



Si— 1 S2 — 1 Sm — '^ 



3 



where 

Eo{si,...,Sm) = EMUHsrr.)-U\s2)U\si)MU{si)U{s2)...U{s^))] (8) 

= EMUism + D/2)...U{S2 + D/2)U{si + D /2))tv{U{s i)U{s 2)... U{s „,))]. 

I. LOWER BOUNDS AND NUMERICAL RESULTS 

In this section we present lower bounds on IA2I based on random walks on a Cayley tree and then provide some 
numerical results. In the Hermitian case, it is possible, for certain choices of Si, ...,Sm in either Eq. ([7|) or Eq. (O, 
that the trace iv[U{si)U{s2)---U{sm)) can be reduced to a trivial trace of the identity matrix by canceling successive 
appearances of U{s)U{s + D/2) and replacing them with 1. The contribution of such choices to Ei is proportional 
to a return probability of a random walk on a Cayley tree as will be seen. 

We begin with an upper bound on the number of such choices: consider the unitaries U{si)U{s2)---U{sk) for some 
k, < k < m. After making all possible cancellations of successive terms, U{s)U{s + D/2), this sequence of unitaries 
may be reduced to another sequence of unitaries U{s[{k))U{s'2{k))...U{s'i^f.-^{k)), for some l{k) < k. Then, consider 
the sequences of unitaries U{si)U{s2)--.U{sk+i)- After making the same cancellations, and then possibly canceling 
U{sk+i) against [/(sj^^,^ (A;)), we find a new sequence of unitaries, U{s[{k + l))U{s'2{k + l)...[/(sj|.^^^^^j (fc + 1)) with 

l{k+l) = l{k)±l and s'j{k + l) = s'^ik) for j < l{k). When l{k + 1) = l{k) - 1, then s^+i is determined by Sk- When 
l{k + 1) = l{k) + 1, then if l{k) > there are D — 1 possible values of Sk+i, while if l{k) — there are D possible 
values. Note that l{k) > for all k. We define N{l{m), m) to be the number of choices of si, S2, s™ which give rise 
to the given l(rn). This is precisely the number of random walks of length m, on a tree with D daughters at the root 
and D — 1 daughters for every other node, that end at a distance l{m) from the root. Note that A^(0,m) is equal to 
D™ times the return probability of a random walk of length m on the Cayley tree. 
An upper bound on A^(0, m) is given by 

^(0, m)<{D- \rl \ 7; <{D~ l)™/22-. (9) 

(m/2)!(m/2)! 

To show Eq. we consider a related problem: consider sequences of HJt) in which may become negative, while 
the number of choices of Sm is considered to be Z? — 1 whenever Z(fc + 1) = Z(fc) + 1, and the number of choices is 
considered to 1 whenever — l(k) — 1. This give an overestimate of the number of sequences, and gives the value 
in Eq. ®. 

On the other hand, a lower bound on iV(0, to) is given by assuming that if l{k + 1) = l{k) + 1 there are only D — 1 
possible choices of Sk+i, regardless of Z(fc), in which case we find that 

N{0,m)>cx{2VD~l)"'/{m + lf/^, (10) 

for some constant c of order unity. These bounds, ()9I10|) . are compeletely standard boundsfi*], and we only repeat 
their derivation for completeness. 

We now use Eq. (jlOp to get a lower bound on IA2I for any completely positive map where the matrices A{s) are 
given by Eq. ([2]). We emphasize that, while the upper bounds elsewhere in this paper are upper bounds on the typical 
behavior of IA2I, the present result is valid for any such map given by Eqs. (|ll2p . and is a quantum analogue of the 

Alon-Boppana boundjlil. First, lAaT = 1 + Ef=2 l-^aT < I + (N^ ~ 1)1^21™ < 1 + iV^|A2r. Note that the 

product of traces in Eq. © is equal to |tr(t/(si)/7(s2)...C/(sm))P and so is positive for all choices of si,...,Sm- If 
l{m) = 0, then \tr{U{si)U{s2)---U{sm))\'^ = N'^, and so the contribution of terms with l{m) = to the sum in Eq. ^ 
is equal to N^N{0,m)/D"'. Therefore, 

a=l 

Thus, IA2I > Ah(c/to3/2)1/'"[i „ ^3/2/(^;^m^2)]i/m^ Picking TO = [ln(ciVV2) - (3/2)ln(ln(ciVV2))]/ln(l/Aij), we 
find 

IA2I > Aff(l-0(ln(ln(Ar))/ln(7V))). (12) 

A very interesting question is to see whether a bound such as (jl2p still holds for arbitrary trace preserving, completely 
positive Hermitian maps £{M). As a partial step towards this more general bound, note that the bound of Eq. (|12p can 
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FIG. 1: Eigenvalues from numerical diagonalization of a completely positive map based on the construction in [3| using expander 
graphs, for A'' — 20, 30, 50. The eigenvalue with eigenvalue unity is not shown. The second largest eigenvalue is at roughly 
\/3/2. Only a single realization is shown for each A''. The inset shows a detail of the behavior at small a. Curves in the inset 
are — 20, 30, 50 from top to bottom; the curves in the main figure are not distinguishable. 



be readily generalized to the following case: let A{s) = ^y P{s)U{s), with the numbers P{s) obeying X]s=i ^i^) ~ 1' 
and with U{s) = U{s + and P{s) = P{s + D/2). Eq. © is a special case of this with P(s) = l/D. 

As stated before, the main result of this paper is that, for any e, when the unitary matrices are chosen randomly, 
the probability that IA2I is within e of \h approaches unity when N becomes large. Interestingly, this seems to be 
true in more generality than just for unitary matrices chosen with the Haar measure. Using the construction in [1], in 
which we pick a random graph with constant coordination number and derive unitary matrices from that graph and 
from certain random phases, numerical studies also show that IA2I is close to \h- We show in Fig. 1 the results of 
numerical diagonalization of systems of size N — 20, 30, 50, so that there are 400, 900, 2500 eigenvalues respectively. 
The second largest eigenvalue is indeed very close to •\/3/2. After sorting the eigenvalues by Aa, from most positive to 
most negative, we plotted the eigenvalues as a function of a/N^: the scaling collapse of the curves is extremely good. 

II. BOUNDS ON EIGENVALUES 

A. Schwinger-Dyson Equations 

We will develop a perturbation theory in 1/A'' to estimate the average (O, which is a product of two traces. To do 
this, we will develop general machinery for computing the average over the unitary group of products of an arbitrary 
number of traces. Consider such a product of the form: 

E[LiL2...Lkl (13) 

where 

Li =tr(C/(si,i)C/(si,2)...C/(si,™J), L2 = tr([/(s2,i)C/(s2,2)...C/(s2,™J), ... (14) 

Here we have an average of k traces, Li, ...Lk, each of which is a product of ruk unitary matrices. We now present 
the Schwinger-Dyson equations. 

Let T", for a ~ l...iV^, be Hermitian matrices such that 

^t;,t;, = <5^.<5,,. (15) 

a=l 
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Then 

Y^iT-T-),, = N5^,. (16) 

a=l 

To compute the average in Eq. (fT^ . we begin with 

S[tr(T°i/(si,i)C/(si,2)...C/(si,„J)L2...ifc]. (17) 

We then use the invariance of the average over the unitary group under an infinitesimal change in variables: 

U{si)^{\ + ieT'')U{si) (18) 
U{si + £1/2) U{si + D/2){1 - zeT"), 

where we recall that U{s + D /2) = U (s)^. Applying the change in variables given in Eq. ([T5)) to Eq. (|17[) . and summing 
over a and dividing by N, we find that: 



S[tr(C/(sia)C/(si,2)...C/(si.„H))i2...ifc] (19) 

, nil 

= -l^E'^^i,-«i..^[t'^(f^(^l^l)-^(^l^^--l))t^-(f^(^l^^-)-C/(si,mJ)i2...ifc] 



i=2 

mi 



^E'5«i^-^i..+-/^^^'"(f^(*i^i)---^(^i.^0)tr(f/(s,+ia)...C/(si,™J)L2...Lfe^ 

i=2 j = l 
^ k mi 



N 

1=2 3 = 1 

We simplify the second and fourth lines after the equality sign of the above equation using U{s)U{s + D /2) = 1 to 
get 

E[iY{U{si,i)U{sia)-U{si^r,^,))L2-Lk] (20) 
1 

= -]^E'^«i.-^i..^NC/(si,i)-C/(si,,-i))tr((7(si,,)...f/(si,„J)L2...Lfe] 
i=2 



-, mi 

+ ]^E^^i.-^i.>-/^^[*^(f^(^i'2)...f/(si,,-i))tr(f/(s,+i,i)...C/(si,„J)L2...Lfe] 

-^^^'^.i,i,s,,,^[tr((7(si,i)...f/(si,„,j;7(.s,,,)C/(s,,,+i)...C/(s,,,_i))i2...ii-i^ 
;=2 0=1 

^ k mi 

J7YYSsi,i,si,j+D/2E[tT{U{si^2)---U{si^„^^)U{si^j + l)U{sij+2)---U{sij^l))L2...Li^lLi+i...Lk]. 



N 

1=2 3 = 1 

These Schwinger-Dyson equations are quite long when written out, but in reality are quite simple. Let us apply 
them to compute the average of tr([/)tr(J7^) over unitary matrices U. We find after one iteration of Eq. ([20]) 
that this is equal to (l/7V)tr(l) = 1. Now consider a more complicated example, to compute the average of 
tr(C/C/)tr([/^[/^) over unitary matrices U. Then, the Schwinger-Dyson equations give after the first iteration: 
(l/iV)(-£:[tr(C/)tr(C/)tr(C/tC/t] + 2£;[tr(C/C/t]) ^ -(l/7V)£;[tr(i7)tr(t/)tr(;7tt/t] + 2. We then apply the equations again 
to the average £;[tr(C/)tr([/)tr(C/tC/t], giving £;[tr(C/)tr([/)tr(C/^C/t] = (l/Ar)(_£:[tr([/[/)tr(C/tc/t)] + 2£;[tr(C/)tr([/t]). 
Since we have already worked out E[tr{U)tr{W]) = 1, we have £:[tr(C/)tr(f7)tr(C/tc/t] ^ ~{l/N)E[tr{UU)tr{WW)] + 
2. Thus, putting it all together, we find that 

E[tr{UU)iT{U^U^)] = 2 + {l/ N^)E[tv{UU)iv{U^U^)] - i2/N^), (21) 
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and hence E[tT{UU)tr{U'fU'<)] = 2. 

We now describe the general algorithm for reducing traces of the form . We initially cancel all pairs of matrices 
U{s)U{s + D/2) appearing successively in the same trace, replacing them with 1. We then apply the equation (|20p . 
Then, we cancel all pairs of matrices U{s)U{s + D/2) appearing successively in the resulting traces, replace the trace 
tr(l) by N, and repeat this procedure for each term. After the first application of Eq. ([20]), the number of terms on 
the right-hand side will be at most rritotai — 1, where mtotai = Tti + m2 + ... + nik- Applying the equations repeatedly 
will generate more and more terms at each application. We regard this as a branching process: each term on the 
right-hand side can be then fed back into the left-hand side of the equation to generate new terms on the right-hand 
side. Note that at every stage, each term will produce at most nitotai — 1 terms on the right-hand side since the total 
number of unitary matrices which appear in the traces, mi + m2 -I- ... -I- m^;, will always be at most mtotai- If the 
number of unitary matrices becomes equal to zero in some term after n iterations of Eq. (PO)) . then we are left with 
only trivial traces and we say that this term "terminates" at level n. 

This algorithm generates an infinite series, where the n-th term in the series is equal to the sum of all terms 
terminating at the n-th level. We claim (and will show later when we discuss the convergence of the series) that, if 
mtotai < N, then this series is absolutely convergent and also the average of the original trace is equal to the sum 
over all levels, n > 1, of the terms which terminate at each level, so that the series converges to the desired answer. 

This series is in fact an infinite series for many simple examples. In fact, for E[tr{UU)ti{U^ U^)] we find that after 
two repetitions of the process, the same average E[tr{UU)tr{W W)] has re-appeared, as can be seen on the right-hand 
side of Eq. (PT|) , and thus the algorithm above does not ever finish because there are always some terms with nontrivial 
traces. In this particular case, however, although the algorithm does not ever finish, the sum of the terms terminating 
at any level n > 1 is equal to zero; in other casesfT^] this is not true and the given series has an infinite number of 
nonvanishing coefficients. We will later see how this infinite series is related to an infinite series in for the given 
trace. 

We win apply this procedure to the trace Eq = E[tT{U{sm+D/2)...U{s2+D/2)U{si+D/2))tT{U{si)U{s2)...U{s„,))]. 
Thus, Li = tr{U{s,n + D/2)...U{s2 + D/2)U{si + D/2)), and L2 = tT{U{si)U{s2).-U{sm))]- Begin by reducing all 
successive pairs of a unitary matrix followed by its Hermitian conjugate. What is left is two traces Li, L2 such that 
TOi = m2 and si^i = S2,m2+i-i + D/2. Wc will proceed by estimating the probability of different values of mi given 
a random choices of si, s„i, and then estimating the behavior of Eq for the given mi = m2. 



B. Length of the Reduced Trace 



In this subsection, we will estimate the number of choices of si, Sm such that the reduced traces, ii, L2 have a 
given mi = 7712. 

We start with the case mi = TO2 = in which case -Eq = N"^ ■ The number of different choices of si, ...,Sm with 
TOi = m2 = is given by Eq. ^ so the contribution of all such choices to Ei is bounded by 

N^D-'^'iD - l)™/22™ ^ N^X'S- (22) 

We can also bound the number of choices of si, ...,Sm which give a given mi > 0. In this case, Z(m) = mi and 
s'j{m) = sij. Using the same argument as gave Eq. the number of such choices is bounded by 

(£)_l)™i/2(£l_l)™/22™. (23) 

This number is independent of the particular values of si.i, si.mi . There are [D/{D — 1)]{D — 1)'"^ different 
possible values of Si_i, Si_mi and therefore the total number of choices of si, Sm which give rise to a given choice 
of si,i, si^rrii is bouudcd by 



D - 1 



(1 \ ^1 

-_) (i?-l)'"/22™. (24) 



D V/D 



C. Nontrivial Words 



We now consider the case the mi > 0. After the first application of Eq. the term on the fourth line with I = 2 
and i = m reduces the trace to (l/7V)i;[tr(f/(si,2)...t/(si,mi)t^(s2,i)---C/(s2,mi-i))] = {l/N)E[l] = 1, so that the 
series terminates at level n = 1. There may also be other terms which reduce the trace to a trival one after a single 
application of Eq. (|20p if the sequence of values Si,iSi,2---Si^mi has a symmetry under a shift: Sj 1 = Sj+m^/o 1 for some 
o > 1 which we refer to as the period of the shift. Here, we treat the index j as periodic with period mi. For example. 



7 



the problem studied in Eq. ((2T|) has such a symmetry under a shift with o = 2. In the event that there is such a shift 
symmetry, then the sum of terms terminating at level n = 1 is equal to o. For a given rrii, the number of choices of 
si^mi which have a shift symmetry with period o is bounded by [D/(_D— !)](£)— 1)™^/°. Thus, from Eq. (p4|) . the 

total number of choices of si, Sm which give rise to a given mi, o is bounded by (D—l)"^^'°(^y==] (_D — 1)™/^2'". 

Thus, the contribution to Ei of terms terminating at level n = 1 is bounded by 



E 



mi 



mi<m 0—2 



o{D-ir^/''[-^2=y\]}. (25) 



The term in Eq. (|25p with o = 2 is is bounded by 2m A^, while the sum of terms in Eq. with o > 2 is bounded 
by {D - l)-i/V[l -{D- 1)"1/S]2A^, which for D > 4 is bounded by 30A^ so that Eq. ([25]) is bounded by 

l + 2mA^ + 30AS. (26) 

We will now bound the sum of all terms terminating at level n > 1. Assuming that the sequence si,2---'Si,mi lacks 
the shift symmetry discussed above, this is the only term which terminates at level 1, and the other terms which 
appear after the first iteration do not terminate and continue to branch, but some of their descendents will terminate 
at lower levels. 

We can estimate the value of a term which terminates at a given level n > 1 as follows. First, there is a sign equal 
to plus or minus 1. Next, there is a factor of (1/A^)". Finally, there is a factor of N for each trace of the form tr(l) 
that appeared in this process. Suppose there are p such traces, giving a factor of N^. How big can p be? Initially we 
have k = 2 different traces. The given term at level n arose from a specific choice of terms on the right-hand side of 
Eq. (PD|) on the first iteration. This specific choice has fci different traces in it, with fci equal to either 1 or 3. After 
the second iteration there are k2 traces, then fca, and so on. The number of traces ^2,^:3, ... can be determined as 
follows: an application of Eq. (|20p may increase the number of traces by one if the term arises from the first or second 
line on the right-hand side, or may decrease the number of traces by one if the term arises from the third or fourth 
line on the right-hand side of Eq. (|20p . Next, some of the traces may be trivial. In the event that the term arose from 
the first, second, or third line of Eq. it is not possible for any of the traces to be trivial, under the assumption 
that any repetitions of the form U{s)U{s + D/2) have been replaced by 1 in the trace. However, in the event that 
the term arose from the fourth line, then it is possible for one of the traces to be trivial, increasing p by one. Thus, 
for each b < n, kb — k^-i is equal to either -1-1, —1, or —2. Let q be equal to the number of times the first or second 
line was used from Eq. (|20p and n — q equal the number of times the third or fourth line was used. Then, in order for 
all traces to be trivial in this particular term resulting from n iterations of Eq. (I20|) , 



2 + q-{n-q)-p = 0. (27) 
Also, since p can only increase when a term from the fourth line is used, 

p<n-q. (28) 

Thus, 

P<l{2 + n)/3\. (29) 
Therefore, the value of a term terminating at the n-th level, n > 0, is bounded in absolute value by 

The number of terms terminating at the n-th level is bounded by 

(2m -1)". (31) 

Note also that there are no terms terminating at level n — 2: if the term does not terminate at level 1, then there are 
either 1 or 3 traces after the first iteration of Eq. pp)) , and then there is no way to have the term terminate at level 
2. Thus, the sum of terms terminating at level n > 1 is bounded in absolute value by 

{2m)^N-^ + {2myN-^ + {2mfN^^ + (2m)^iV"'* + (2m)^A^~^ -t- ... (32) 

< Sm^iV"^ + lem-^A^^^l" i _^ 2mN-^ + Am^N^A \ ^. 

V /I — 6riv^N~^ 
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D. Convergence of Series 

We now show the claim that, for rritotai = + ^2 + •■• + '^fe < N, the average E[LiL2...Lk] is indeed equal 
to the sum over all levels n > 1 of the number of terms terminating at each level and that the series is absolutely 
convergent. After n iterations of Eq. (|20p some of the terms have terminated. There are at most {rritotai — 1)" terms 
which have not terminated, since there are at most {rritotai — 1)" terms. Each of these terms is equal to plus or minus 
one times A^^" times N^" where p„ is the number of times a trivial trace appeared in the process, times the average 
of a product of traces. There are at most rritotai — Pn different traces in the product since there were originally at 
most rritotai unitary matrices. Thus, since each trace is bounded in absolute magnitude by N, the sum of all terms 
which have not terminated after n applications of Eq. ([20]) is bounded in absolute value by {rritotai — 1)"7V~"A^™*°'''' 
which converges to zero as n — > cxi for mtotai < N . Thus, the difference between the sum of the terms terminating 
at the first n levels and the actual value of the average E[LiL2...Lk\ converges to zero as n — > 00. The sum of all 
terms terminating at a given level is bounded in absolute value by the number of such terms, times N~"'NP" , and so 
is bounded by {rritotai — l)"iV~"A'''"*°*"' and so the series is absolutely convergent for rritotai < N. This shows the 
desired claim. 



E. Loose Bound 



Adding the results in Eq. (|22I26I32|) . we find that for 2m < A^, 

El < 1 + [N^ + 2m + 3Q]X'Ji + 8m^N-^ + 16m'^N-^(l + 2nrN-'^ +4m^N-A Krr^- (33) 

V /I — 8m'' A* ~^ 

We now pick m = log(A^'*)/log(l/A^f), so 

El <l + 16{l + o{l))[\og{N^)/\og{l/XH)]^N-^, (34) 

where o(l) denotes terms asymptotically tending to zero as A^ ^ 00. Thus, the average of IA2I over the unitary group 
is bounded by 

r log(l/Aff )/log(A''') 

{l6(l + o(l))[log(A^4)/i^g(^/^^)]4^-2| (35) 

= {l+o{l))\ioose{D), 

where 



Xioose{D) = = y . (36) 

Further, using Markov's inequality, the probability that IA2I is greater than cXioose{D)j for any c > 1, is bounded 
by (1 + o(l))c~'°s(^*)/'°s(i/^ff)^ so that for large A^ it is very rare for IA2I to be significantly above the loose bound 

Mooae{D). 



F. Tight Bound 



We now tighten the bound. On a given iteration of the Schwinger-Dyson equations, we go from a product of k 
traces to a product offc + l,fc— 1, orfc — 2 traces. We will keep track of how the matrices move under this iteration 
process using a function /«((') *)) from pairs of integers to pairs of integers. We say that the matrix U{si^i) in the 
given product of traces, LiL2...Lk, is in position {l,i). Let us consider the case of a term on the first line, where m 
increases by one. Then, for any given j in the sum on the first line, we say that the matrix in position (i, 1), for i < j 
on the n + 1-st iteration corresponds to the matrix in position (l,i) on the n-th iteration, and so /„((!,«)) = (1, i), 
while the matrix in position (2, i) on the n + 1-st iteration corresponds to the matrix in position {l,i + j — I) on the 
n-th iteration, so /^((l, i + j — I)) = (2, i). The matrix in position (/, i), for 2 < ^ < fc + 1 on the n + 1-st iteration 
corresponds to the matrix (/ — 1, i) on the n-th iteration, so /„(/ — 1, i) = i). We follow a similar procedure for the 
other lines of Eq. (PO)) and if there are cancellations, we keep track of how the matrix moves under the cancellations. 

We then keep track of which matrix after n iterations corresponds to a given matrix before any iterations, by defining 
Fn{{l, i)) = /„(/„_i(.../i((/, i)) for I = 1,2. Let us say that the matrix at position {I, i) for Z = 1, 2 is "trivially moved" 
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under the n-th iteration of the Schwinger-Dyson equations if we are considering a term in the equations which did 
not arise from T°[/(s; i)); that is, the matrix is trivially moved if it is not in position using a term on the 

first or second line, or in position using a position from the third or fourth line, or in position (1, 1). Let us 
define a "rung cancellation of matrix i" to be the case in which, for some n, after the n-th iteration of the Schwinger- 
Dyson equation we perform a series of cancellations such that the following hold[13]. First, a matrix in position 
is canceled against a matrix in position {l',j') such that = F„_i((l,i)) and = Fn-i{{2,mi + 1 — i)). 

Second, at all previous iterations up to the n — 1-th iteration, the given matrix was trivially moved. If there is a 
rung cancellation of matrix 1 on the first iteration, then all matrices cancel and the trace is equal to unity; this is 
precisely the case with I — 2, j — m discussed at the start of the section "Nontrivial Words". Note that the matrix in 
position — i^„_i((2,TOi -|- 1 — i)) is equal to C/(s2,mi+i-i) — U{si + D/2) which is why the matrix in position 

{l,j) — F„_i((l,i)) can be canceled against this matrix. 

We now make a stronger claim: for any given i, the sum of all terms with a rung cancellation of matrix i 
is equal to unity. To show this, consider the trace tr(C/(Sm + D/2)...U{si+i + D/2)X^U{si-i + D/2)...U{si + 
D/2))tT{U{si)...U{si-i)XU{si-\-i)...U{sm), where X is some arbitrary unitary matrix. Averaging this trace over 
all unitary matrices U{s) and over all unitary matrices X with the Haar measure, we find that the trace is equal to 
unity. However, applying the Schwinger-Dyson equations to this trace generates precisely the sum of terms mentioned 
above, those in which there is a rung cancellation of matrix i. Thus, this sum equals unity. We further claim that for 
any given ii, ^2, id, the sum of all terms with rung cancellations of matrices ii, 12, ■■■id is equal to unity, as may be 
shown by considering a trace in which matrices [/(s^j), U{si^), ... are replaced by Xi,X2, and the trace is averaged 
over the different Xi, X2, .... 

On the other hand, if a term terminates at level n and matrix i does not have a rung cancellation, then at some 
previous iteration n either the matrix was not trivially moved or was canceled against a matrix in position j) such 
that = Fn-i{l' with ^ (2, mi -1- 1 — i). In the later case, for I' = 2 we know that Sm+i~j' + D/2 = Si, 

while for I' = 1 we know that Sji = Si, thus in both cases identifying some k ^ i such either si^i = si^k or 
si.i — Si, A; + D/2. If the matrix was not trivially moved, we can also identify some k ^ i with the same properties. 
Let us write k = T{i) in both cases, for some function T(j). 

Now consider the sum of terms in which for no i is there a rung cancellation of matrix i. By the inclusion-exclusion 
principle in combinatorics, this is equal to the sum of all terms, minus the sum over i of the sum of terms in which 
there is a rung cancellation of matrix i, plus one half the sum over ii ^ 12 of the sum of terms in which there are 
rung cancellations of matrices ii, 12, and so on. This is equal to the sum of all terms minus the sum 

mi — mi(mi — l)/2 + mi(mi — l)(mi — 2)/6 — ... = 1. (37) 

Thus, the sum of all terms is equal to one plus the sum of terms in which for no i is there a rung cancellation of 
matrix i. So, we now focus on the sum of terms with no rung cancellation, which we define to be Eq[si, Sm)- If 
si, Sm has a shift symmetry as above, then there may be terms in this sum terminating at the first level; the sum 
of these terms is o — 1. 

Each Ef) we are averaging over the unitary group results from a particular set of choices of si, in the sum 
in Eq. ([T]). There are different terms in this sum in ([7]). We begin by bounding, for any given level n, the 
number of choices of si,...,Sm which give rise to an Eq which produces a term in the Schwinger-Dyson equations 
which terminates at level n with no rung cancellations. Suppose for a given choice of si, Sm there is such a term 
which terminates at level n with no rung cancellations. There were two traces of m unitaries in the definition of 
Eq] then, after canceling successive pairs U{s)U{s + D/2), we have mi < m unitaries for some mi. The number of 
different initial choices of si, S£) which produce a given mi after these cancellations is given in Eq. (j24p . 

Then we iterate the Schwinger-Dyson equations with a particular choice of l,j at each level, where 1 <l < 2mi and 
1 < j < 2mi as given in Eq. if we pick a term from the first or second line of the Schwinger-Dyson equations, 

we set I = \ aX that level. At each such iteration of the Schwinger-Dyson equations, there may be cancellations in 
two different traces if the term came from the second line of Eq. (j20p . with at most mi cancellations in each trace, or 
cancellations in two different places of a single trace, if the term came from the fourth line of Eq. (|20p . with at most toi 
cancellations in each place. Let us call the number of cancellations ci, C2 with < ci < mi and < C2 < mi. Then, 
by specifying ci, C2 for each iteration, we succeed in fully specifying how the matrices move under the n iterations 
of the Schwinger-Dyson equation; this requires specifying 2n numbers ranging from 1...2mi, and 2n numbers ranging 
from 0...mi. In particular, since there are no rung cancellations, we succeed in specifying for each i, 1 < i < mi, some 
i ^ i such that either si^i = Sij or sis — si,j + D/2, giving the function T(i). Having specified this function, there 
are now only at most [D/{D — !)](£' — i)'"i/2 possible values of si^i, si^m^. To show this, we start by specifying 
the value of si_i, which can assume any of D different values. By specifying si^i we have fixed the value of of Si^^(i), 
as well as the value of any j such that t(j) = 1, so that there are now at most mi — 2 different values of si^i which 
remain undetermined. We then find the smallest ji such that sij-^ is undetermined and specify its value. Note that 
there are only D — 1 possible values of this sij^ since, by assumption, sij^ ^ sijj_i -I- £)/2. Having specified this 
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sijj, we have fixed the value of Si^r(ji) as well as the value of any j such that T{j) = ji. We then find the smallest j2 
such that si j2 is undetermined and specify that value. Proceeding in this way, we succeed in specifying si i, Si^mi 
by specifying one of at most [D/{D — 1)]{D — l)™i/2 different choices. Thus, there are at most 

[D/{D - !)](£> - l)™i/2(2mi)2"(TOi + 1)2" < [D/{D - 1)](L> - l)"''^\2mi)'^" (38) 

such choices of si,...,s„ij which can produce a term which terminates at level n. Using Eq. (j24[) . the number of 
choices of Si, which can produce a term which terminates at level n is at most 

^ /Oil \4n+l 

V (D - l)'"/22'"(2mi)4" <(D- l)"'/^2" ^ ' . (39) 

An + I 

mi— 

For any Si,...,Sm, we define ?^mm(si, Sm) to be the smallest level at which a term terminates with no rung 
cancellations. We re-write the sum in Eq. ([7]) as 

^ ^ oo D D D 
^l = ^+(n)"E E E - E ^n^..isu...,^^hnKi^l,-,Srn), (40) 



n— Osi — 16*2 — 1 Sti 



so that the second sum is over the set of all values of Si, with the given rimm = We note that the bound of 
Eq. (|30p continues to apply to the terms terminating with no rung cancellations, and the bound of Eq. pip continues 
to bound the number of such terms terminating with no rung cancellations. From Eq. (|30p. a bound on the value of 
the term terminating at the n-th level, for any n > is 

Therefore, for any si, s^, 

i?^(si,...,s™) < iV^^'/'^"(2m-l)" (42) 

^ ^^,2 [iV-^/^(2m-l)]""^- 
1 - Ar-2/3(2m- 1) 

From Eqs. (|39|40|42p . 

p ^ 1 , .,2^,nv (2m+l)^"+i [N-y^{2m-l)r 

E, < 1 + N X^l^ ^^^^ l-N-V^{2m-l) ^^^^ 



oo 



■^^^ (4n+l)[l-7V-2/3(2TO-l)]^ 
We then pick m = (l/4)7V2/i5^ go th^t A^-2/3(-2m + 1)^ < 1/2 and 

IA2I < (£;i-l)i/'"<iV'/™AH(l + 0(l))i/'" (44) 
^ \H{l + 0{\og{N)N-^/^^). 

As before, using Markov's inequality, the probability that IA2I is greater than cXh{D), for any c > 1, is bounded by 

C-(l/4)W=/^=(l + 0(l0g(7V)iV-2/15). 

This shows that for any e, the probability that A2 < Xh + e approaches unity as — > cx). Combined with the 
previous lower bound (|12p . this proves the main result. 



III. DISCUSSION 



We consider some analogies between these results and lattice gauge theory, some applications of these results, and 
some extensions. We begin with analogies between the random construction of quantum expanders and lattice gauge 
theory and the Eguchi-Kawai construction p^. 
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A. Gauge Theory Analogies 

Consider a lattice gauge U{N) theory in D/2 dimensions on a hypercubic lattice, with unitary matrices U^{x) 
defined for each link of the lattice. Here, x represent a point on the lattice, and d represents the direction of the link: 
if d < D/2, it points in the direction of increasing the d-th coordinate by unity, while if D/2 < d < D, then it points 
in the direction of decreasing d — D/2-th coordinate by unity. Then, for a given choice of si, Sm we can define a 
path, starting at the origin, and then moving in direction Si,S2, ■■■ until rn steps have been taken. We can define a 
product of traces associated with this path: tr([/sj (0)J7s2 (0 + si)...)tr(...C/s2 (0 + si)'''t/si (0)'''). For certain choices of 
the si, ... this path returns to the origin after m steps, in which case the product of traces is a product of two Wilson 
loops. If, however, the path does not return to the origin, the product of traces is not invariant under non-Abelian 
gauge transformations, and hence the average of the product of traces is equal to unity. 

At infinite coupling, all of the unitary matrices are independent, except for the constraint U^{x) = U^^jy^^{x), and 
even if the path does return to the origin, the average of this product of traces is equal to unity, unless, by chance, the 
path of length m exactly retraces itself. The probability of this retracing, for a random path, is precisely the Cayley 
tree return probability discussed previously. Thus, this lattice gauge theory at infinite couphng has tr([/sj (0)t/s2 (0 + 
Si)...)tT{...Us2 (0 + si)^Usi (0)^) = 1 + N^D~^^N{0, m). The Eguchi-Kawai construction is an approximation to large 
N gauge theory which replaces the infinite lattice by a single site: this turns tr([/si (0)L/s2 (0 + si)...)tr(...C/s2 (0 + 
si)^C/si (0)^) into £^o(si, S2, ...), the quantity we considered before. Thus, this paper can be seen as an estimation of 
corrections to the Eguchi-Kawai construction in the infinite coupling limit. There are a number interesting terms in 
these corrections: for example, the average tr(C/(l)C/(l))tr(J7^(l)C/^(l)) is equal to 2 as calculated before, but the 
corresponding average in the lattice gauge theory is equal to 1. 



B. Applications 

The general properties of ground states of local Hamiltonians with an excitation gap have become of great interest 
recently. A basic result [l^ [isj is that correlations decay exponentially in such systems. One application of quantum 
expanders is to finding matrix product states of one-dimensional quantum systems with the following properties: the 
correlation length is of order unity, the Hilbert space dimension on a single site is small, also of order unity, and yet 
the entanglement entropy across any cut is large. As an example, consider a matrix product state of the form: 

*(si,S2,---SAr) = ^ Aq,/3(si)A^,^(s2)A,5(s3)--- (45) 
a, 13,... 

where si, S2, sat are spin variables in a one-dimensional quantum system of sites. Associated with the matrix 
product state is a completely positive map as in Eq. ([1]). If this map has a gap in its eigenspectrum to the second 
largest eigenvalue, then the state \1/ has exponentially decaying correlations [ig, so that if operator A has support on 
sites 1, j and operator B has support on sites j + 1, N , then (^, A_B5') — (^I*, B'^) is exponentially small 

in as required for the ground state of a gapped, local quantum system. However, as discussed in pj], this means 
that the existence of quantum expanders implies that the mere fact of exponentially decaying correlations does not 
suffice to prove bounds on entanglement entropy. Instead, bounds on the entanglement entropy[l3| proceed through 
a different route and currently give weak bounds. 

However, in (T7| . a conjecture was developed regarding properties of quantum expanders that may help in prov- 
ing tighter bounds on entanglement entropy. Consider the following different correlation function. Let A be an 
operator with support on the sites 1,2..., j — I and j + l,j + I + 1,...,N. Let v]/ = ^^^^ A{a)'i' L{a) (8) ^'_R(a), 
where ^'^(a) are orthonormal states on sites and ^_R(a) are orthonormal states on j + 1,...,N. Let 

Bl = '^('^)^i,o(ci^))(^-L,o <8i 1_R, where 1_r, is the unit operator on Xj+i^n, for some function 0{a). Then, 

it was shown that for a gapped local Hamiltonian, 

(*o,^5l*o> - (*o,A*o>(^'o,Bl*o) < \\A\\\\B\\0{exp[-l/lo]), (46) 

for some Iq. 

It was conjectured in [13] that there is a function f{Deff) such that if Eq. (|46p holds for a state ^E" for some Iq, 
then the entanglement entropy of ^E" is bounded by f{D^°). Interestingly, it seems that an expander where the A{s) 
are random unitaries is unlikely to satisfy Eq. (|46p . If this could be shown to be a general property of expanders, 
showing the conjecture, this would provide another way of studying area laws in quantum systems. 
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C. Extensions 



The method of Schwinger-Dyson equations used here is fairly general and could be applied to other groups, such 
as 0{N) or Sp{2N). We have not done the calculation, but it seems that random choices from these groups will also 
give quantum expanders. Always, the unit matrix is an eigenvector of the map £{M) with eigenvalue unity. Any 
matrix in the center of the group is also an eigenvector of £{M) with eigenvalue unity, but for these cases, all elements 
of the center are proportional to the identity matrix, and thus do not give rise to additional eigenvectors with unit 
eigenvalue. 

The method can be directly extended to the non-Hcrmitian case. Some of the combinatorics become slightly easier 
here. From Eq. (|3]), the average of I'^aP™ O'^^r the unitary group is bounded by the average of the trace: 

(^)'"Ef,=i-Ef„=iE^,=i-El=ii?[ tr{U{-s,)U{-S2)...U{-Srn)U\s„,)...U\s2)UHsi)) (47) 

y^ti{U{s^)U{s2)...U{s^)U\-Sm)-U\-32)U\si))] 

The probability of having U {si)U {s2)...U {sm)U^ {s,n)...U^ {s2)U^ ^i) cancel to the identity matrix is equal to 1/D™. 
Let 

A„« ^ ^. (48) 

Carrying through the calculation one finds that, for any e > 0, the probability that IA2I < XnH + e approaches unity 
as iV — > 00. Note that XnH < Xh and also note that in the non-Hermitian case even a tight estimate on the average 
of the trace only provides an upper bound on the eigenvalue, due to the inequality in Eq. ([4]). However, numerical 
work suggests that the eigenvalue is asymptotically equal to X^h with high probability in this case. 

We can also provide a lower bound on the trace in the non-Hcrmitian case. For any choice of unitaries U{s), the 
sum of terms in (|T7)) with s,; = Si is bounded below by 1/D™, while the other terms are all positive, so that the given 
average (|47)) is bounded below by N^D~"^. This result extends readily to an arbitrary choice of A{s), constrained 
only by the trace-preserving condition, Es=i A{s)A^{s) = 1. 

Acknowledgments — I thank T. J. Osborne for many useful discussions on this topic. This work was inspired 
by discussions with A. Harrow, T. J. Osborne, and F. Verstraete at the workshop on "Lieb-Robinson Bounds and 
Applications" at the Erwin Schrodinger Institute. This work supported by U. S. DOE Contract No. DE-AC52- 
06NA25396. 



APPENDIX A: QUANTUM EDGE EXPANDERS 



In the appendix, we discuss the relationship between quantum expanders and another concept, a quantum version 
of an edge expander. We define a map to be a "quantum edge expander" if the following condition holds: for any 
A^-by-A^ Hermitian matrix P such that — P and such that P has I non-zero eigenvalues, I < N/2, 

triP£{P)) < Aetr(P), (Al) 

for some Ae less than one. We then prove a relation between Ag and | A2 1 , showing that a quantum edge expander is a 
quantum expander. This is a quantum analogue of a theorem of Tanner ,18J and Alon and Milman^lQ], which shows 
that an edge expander has a spectral gap. We consider only the Hermitian case in this appendix, leaving the behavior 
in the non-Hermitian case open. We also assume that the second largest eigenvalue is positive; the case where it is 
negative can be considered by looking at the square of the map £(M). 

Let X be the eigenvector of £ with eigenvalue A2. We work in a basis in which X is diagonal. 




and such that ei > 62 > ... > e^v. Since X is orthogonal to the unit matrix, using the inner product (A, N) — tr{XN), 
we have tr(A) = 0. Define m such that ei > for i < m and Ci < for i > m. Without loss of generality we may 
suppose that m < N/2, as otherwise we could have considered the matrix —X which has the same eigenvalue. Define 
M{i,j) to be the matrix with an a unit entry in the i-th row and j-th column and zero everywhere else. Define 



P,,=tr{M{i,i)£{M{j,]))). 



(A3) 
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Then, since the map £ is trace preserving, we have 

i 

for all j. Also, we have = Pji. Finally, since £ is completely positive, we have > for all z, j. Then, 

N N 

A2 = tr(X£(X))-^^e,e,P.,. 

1=1 j=i 

Define by 



^ e2 i < m 
fi = i> m. 



(A4) 



(A5) 



(A6) 



Then, 



Then, 



i=i i=i 



i— 1 j>i i— 1 j>2 



< 



< 



(l/2)EE^^^2(/2 + /|) 



= V2 



\ 



E E (/^ - /^)' V2v/(r^. 

i— 1 j>z 



(A7) 



(A8) 



where the first inequality uses Cauchy-Schwarz, the last equality uses Eq. (|A4p . and the last inequality uses Eq. (|A7p 
Let Pi be the projector onto the vector space spanned by the first i eigenvectors of M. Then, 



i=l jyi i=2 



(A9) 



4=1 



Using the property of a quantum edge expander, (|Aip , we have 

m mm 

m 

= E/'(l-^e) = (l-Ae). 

Combining Eqs. (|A8IA9IA10p . we find that 



i=l 



1 - A, < V2(l- A2). 



(All) 



We finally show the converse result, that a quantum expander is a quantum edge expander. The normalized 
eigenvector with unit eigenvalue is vi = (l/ViV)l- We have {vi,P) = tr(P)/V7V- So, P = tT{P)l/N + P' , where 
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P' = P-tr(P)l/7V. Thentr(P£:(P)) < \\2\tr{P'^)+{vi, P)^ = |A2|(tr(P)-tr(P)VA^)+tr(P) VA^. Iftr(P) = I < N/2, 
then tr(P)2 < (Ar/2)tr(P) and tr(P£:(P)) < |A2|tr(P) + (1 - |A2|)tr(P)/2 so 

Ae < IA2I/2 + I/2. (A12) 
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